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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3513 


HEAT TRANSFER AT THE FORWARD STAGNATION POINT OF BLUNT BODIES 
By Eli Reshotko and Clarence B. Cohen 


SUMMARY 

Relations are presented for the calculation of heat transfer at 
fonvard stagnation point of hoth two-dimensional and axially sym- 
metric blunt bodies. The relations for the heat transfer^ which were 
obtained from exact solutions to the equations of the laminar boundary 
layer, are presented in terms of the local velocity gradient at the 
stagnation point . These exact solutions include effects of variation 
of fluid properties, Prandtl number, and transpiration cooling. Exam- 
ples illustrating the calculation procedure are also included. 


S INTRODUCTION 

o . 

In the design of supersonic air vehicles, probably the most crit- 
ical areas with regard to aerodynamic heating will be in the vicinity 
of the forward stagnation points of bl\int bodies. This problem can be 
treated using laminar-boundary-layer theory. Stagnation-point heat- 
transfer relations are presented herein based on exact solutions to 
the laminar-boxmdary-layer equations with external flows of the Falkner- 

Skan type (ug = Ax^) and with constant surface temperature. The latter 

restriction is required because the existing exact solutions that in- 
dude surface temperature variation apply to a temperature variation of 

the form (t^ - tgp) = Bx^ . (All symbols ane defined in appendix A.) 

Using this variation it can be shown that the available exact solutions 
for variable surface temperature are not valid at a stagnation point be- 
cause the term ^ ^k has been omitted from the energy equation. 

The applicable solutions may be divided roughly according to their 
consideration of fluid properties. 

Constant fluid properties . - For the case of constant fluid proper- 
’ll® s>SquiiT7rdri)lin^ (ref. 2) have evaluated the parameter 

Nu/-VRi^ Tot l3Vo-dii2i6risioria.l and. axially synmiG'tric ilow^ rGspGctivGly . 
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These calculations, for Prandtl numbers from 0.6 to 1.0, indicate that 
the following approximate relation closely represents the Prandtl num- 
ber effect on Nu/^Rey-: 



Schuh (ref. 3) performed similar calculations with Pr = 0.7 for the 
flat plate, the two-dimensional stagnation point, and an intermediate 
pressure gradient. Levy (ref. 4) calculated many cases of constant- 
property flows for Prandtl mmibers 0.7, 1.0, 5, and 10. Solutions for 
porous surfaces were obtained in reference 5 to estimate the effects of 
transpiration cooling on heat transfer. A numerical comparison of some 
of the results of references 3 to 5 is given in appendix B of reference 5. 

Variable fluid properties . - References 6 and 7 present solutions 
for Pry =0.7 with and without transpiration cooling. The fluid prop- 
erties were assumed to vary as powers of the temperature so that the 
Prandtl number varied within the boimdary layer. Solutions to the com- 
pressible laminar-boundary-layer equations for Pr of 0.7 and 1.0 with 
no transpiration cooling are presented in reference 8. The viscosity 
was taken to vary linearly with temperature so that pp = constant . So- 
lutions are presented in reference 9 for Pr = 1.0 with no transpiration 
cooling. The viscosity law there assumed is similar to that of Chapman 
and Rube sin (ref. lO). 

For the present paper, additional stagnation-point solutions (both 
two-dimensional and axially symmetric ) have been obtained for Pr = 0.7 
with and without transpiration cooling. The method by which these solu- 
tions were calculated is indicated in appendix B. 

The heat-transfer relations are given in terms of (l) a heat- 
transfer parameter Nu/-^/Rey from the exact solutions and (2) the 
local velocity gradient at the stagnation point. 


HEAT-TRAHSFER RELATIONS 


The characteristic quantity representing heat transfer which is 
obtained from exact solutions of the laminar-bo^lndary-layer eqviations 


is Nu/-^Rey, where the Nusselt number Nu = 


'w hx j . V, 

T TT- = and the 

^aw " ^w ^ 


Rey = 


UgX 


w 


Reynolds number 
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In describing the flow over blunt bodies, the first term of the 
Taylor series for surface velocity is the linear relation 

Ug = Cx (2) 


where C is a constant representing the local velocity gradient exter- 
nal to the boundary layer at the stagnation point. Because of this lin- 
ear velocity relation (eq. (2)), the ratio Nu//\/Re^ at the stagnation 
point (x = 0) is finite, althovigh the Nusselt and Reynolds nimbers indi- 
vidually become zero at x = 0. Furthermore, for a stagnation point, 
the external Mach number is zero, so that in the definition of Nu the 
adiabatic wall temperature is the free -stream stagnation tempera tirre tQ 

From the preceding relations, the heat-transfer coefficient h can 
be expressed 



(3) 


Similarly, the local rate of heat transfer q to the wall per unit wall 
area, defined by the relation q = h(tQ - t^), is 


Q. = 


K<.^o - 

/ Nu \ 


J 


(4) 


To solve a given problem, it is necessaiy to determine the velocity gra- 
dient C from potential flow relations, and the parameter Nu/^^Re^ 
from the pertinent bo\mdary- layer solutions . 


Local Velocity Gradient 

The local velocity gradient in the vicinity of the stagnation point 
should wherever possible be evaluated from experimental velocity or 
pressure distributions for the configuration being considered. However, 
in the absence of such e3q>erimental information, a reasonable estimate 
may be made from potential flow relations . 

The velocity on the surface of spheres and cylinders in subsonic 
potential flow is determined in references 11 and 12, respectively. 

The method used therein is the Rayleigh-Janzen expansion in powers of 
Mach number. The resulting velocity gradients at the stagnation point 
f for y = 1.4 are: 

Sphere : 

C = — (1 - 0.252 M£ - 0.0175 ^^) 


( 5 ) 
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Cylinder : 


2u-| p 

C = — (l - 0.416 mJ 


- 0.164 M^) 


( 6 ) 


where and M-j_ are the subsonic velocity and Mach number ahead of 

the body, and r is the radius of the body. Relations (5) and (6) were 
derived for uniform subsonic flow approaching the bodies. For the case 
of supersonic flight speeds, the body is preceded by a curved bow wave 
behind which the flow is nonuniform. Nevertheless, for a sphere the ex- 
perimental pressure distribution near the stagnation point has been 
fovind from reference 13 to be well represented without consideration of 
shock curvature. In this case, up and Mp are taken to be the values 
of velocity and Mach number just behind the normal portion of the shock. 
However, for a cylinder in supersonic flight with axis normal to the 
flow, eqiiation (6) does, not appear adequate. The experimentally ob- 
served velocity gradients (refs. 14 and 15) more closely resemble that 

ui 

predicted for Newtonian flow, which is C = . 


Charts of Heat -Transfer Parameter 

Values of the heat-transfer parameter Nu/^Rev for forced- 
convection heat transfer are presented in figure 1. The solid lines 
are for Pr of 1.0 and 0.7 and represent the solutions of reference 
9 and the present report. The dashed lines are for constant Prandtl 
nxmnber determined using eq^uation (l). The small difference between 
the solutions of reference 7 and those of the present report for 
Pr = 0.7 arises from the different assumption of fluid-property var- 
iation. These solutions show that, for t.^^/tQ between 0 and 2, the 

variation in Nu/ YRS w less than about 10 percent from the values 
for t.^^Ao = 1 for both the two-dimensional and the axially symmetric 
stagnation points. The values for t^/tQ = 1 are also those obtained 
for constant fluid properties. 

The effect of transpiration through a porous wall on the heat- 
transfer parameter is shown in figure 2. The abscissa is the dimen- 
sionless normal velocity parameter which stems from the ex- 

act solutions. In general, as the normal velocity into the boundary 
layer is increased, the heat-transfer parameter decreases significantly. 
This decrease becomes more pronounced as the surface temperature is low- 
ered. The broken lines represent values obtained by interpolation from 
the solutions of reference 7. The differences between the solutions of 
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reference 7 and those of the present report^ as evidenced hy the lack of 
imiformity in the spacing of the curves of figure 2, are due to the dif- 
ferent assumptions of fluid-property variation. For constant fluid 


properties 


S ■ ‘) 


the two calculations shoiild agree. An estimate of 


the spacing between the curves of t^/tQ = 0 and 0.25 for the axially 
symmetric case may be obta ined by examination of the presently calcu- 
lated points for v^/ /\/v ^ = 0.707. 


APPLICATION 

To calculate either the heat-transfer coefficient h from equa- 
tion (3) or the local heat-transfer rate q at a stagnation point from 
equation (4), it is necessary to know the local velocity gradient C at 
the stagnation point and the dimensionless heat-transfer parameter 

Nu//^ew The local velocity gradient may be obtained from established 
experimental data , and Nu/ f\[R^ chosen from figures 1 or 2 according 
to the conditions at the stagnation point. In the absence of experimen- 
tal data, the local velocity gradient may be estimated from equation (5) 
or (s), whichever is appropriate. 

The problem at hand may often, however, be that of finding the sur- 
face temperature \mder the action of various types of cooling or else 
that of finding the transpiration mass flow required to maintain a spec- 
ified siarface temperature. For such problems it is necessary to con- 
sider the usixal heat balance 

^ Pw^w(^)tc + 1r (7) 

The left side of equation (7) represents the rate of heat transfer by 
convection to the sirrface (eq. (4)). The first term on the right side 
qj^^, which represents the rate of heat removal from the surface by an 

internal coolant, is the product of the internal coolant mass flow per 
unit cooled surface area and the enthalpy rise of this internal coolant 
per xmit mass. The second term on the right side of equation (7) ac- 
counts for the heat removed by transpiration cooling, where (AH)tc is 
the enthalpy rise of the transpiration coolant per \mlt mass. The last 
term of equation (?) represents the heat transferred from the surface 
by radiation. This may be quite important when high surface tempera- 
tures are being considered. 
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The solutions shown in figure 2 are exact only if the coolant and 
the external fluid at the stagnation point have the same physical prop- 
erties . In some applications it may be desirable to use a coolant which 
is a different fluid from that flowing over the body {e.g., water as a 
coolant, with air flowing over the body) . Although it is not possible 
from the present calculation to predict the additional effects if dif- 
ferent fluids are involved, it is expected that the present results 
would be reasonable so long as either the rate of injection is small, 
or the properties of the coolant are similar to those of the external 
flow. 

For the general case of a liquid coolant being evaporated. 


The subscripts v and I refer to the vapor and liquid states, respec- 
tively, Ly is the latent heat of vaporization of the coolant, ty is 

the vaporization temperature of the coolant, and t^es is the tempera- 
ture of the coolant reservoir. 


retained and the calciilation performed to effect a heat balance. The 
following form of equation (v) is convenient: 


Several examples are presented to illustrate the use of the heat- 
transfer relations and the heat balance. The calculations are for a 
hemispherical-nosed body having a 1-foot radius and moving at M = 5 
at an altitude of 175,000 feet where the ambient temperature is 170° F. 
The Mach number and velocity behind the shock are Mj_ = 0.415 and 
= 1231 feet per second. Equation (5) indicates a velocity gradient 
at the stagnation point of magnitude 1767 per second. 

I. Rate of internal cooling required to maintain a given wall tem- 
perature. - Since o nly the heat-transfer rate is required, the calcula- 
tion for example I can be performed using equation (i) and the values 
of the heat-transfer parameter Nu/^Re^ fi’om figiire 1. A value of 
q is obtained for each assimied wall temperature. The result is shown 
on the left in figvire 3. 


(^)tc — °p,v^^w ” ^v^ + + ^p,Z^^v “ ^res^ 


For a specified problem, the pertinent terms of equation (7) are 



( 8 ) 
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II. Rate of transpiration cooling required to maintain a given wall 
temperature . - For example II, the heat balance from equation (8) is 


Nu 

‘ P^w(^)tc 

^w 


^p,w("^0 ” ^w) 


For given values of free-stream stagnation temperature and coolant res- 


ervoir conditions^ the bracketed term is a function of wall temperatiJire 
only. This equation thus provides a relation between the heat-transfer 
parameter Nu//^Re^ and the normal velocity parameter v-^^/^v ^ . An- 
other relation between these two q\iantities is that of figure 2. The 
resialts of simiiltaneously satisfying both relations are shown in figure 
3. The coolant assumed was water at a reservoir temperature of 40° F. 

III. Surface equilibrium temperature considering only the heat 
transfer by radiation . - For example III, the heat balance, from 
equation (7), is 



where a is the Stefan-Boltzmann coefficient and e is the emissivity 
of the wall material. The radiation receiver temperattire to the foiarth 
power was neglected here, as it woiild be small compared with . The 

result of the indicated calculation is that, for an assumed emissivity 
of 0.6, radiation lowers the wall temperature from a stagnation temper- 
ature of 3780^ R to a temperature of 1740^ R. 

IV. Rate of transpiration cooling required to maintain a given wall 
temperature with radiation . - The heat balance for example IV can be 
written, from equation (8), 



Here again the coolant was taken to be water at a reservoir temperature 
of 40^ F, and the emissivity was taken to be 0.6. For each assumed val- 
ue of surface temperatiire (necessarily below 1740^ R because of the ra- 
diation calculation of example III), values of Nu/^Re^ and v^/n^v^C 

are found by trial-and-error solution using figure 2. The results are 
shown in figure 3. 
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Discussion of examples . - An examination of figure 3 shovs that, 
omitting radiation, for small reduction in surface temperature the mode 
of cooling does not significantly affect the heat-transfer rate. How- 
ever, as the desired surface temperature is further lowered, transpira- 
tion cooling requires noticeably less heat removal than does internal 
cooling. This same trend has been noted for zero pressiore gradient 
(ref. 16). On a weight-of-coolant basis, transpiration cooling using 
a liquid coolant which becomes vaporized has a significant advantage 
over internal cooling for which boiling may not be permitted. 

Radiation alone is seen to cool the surface from 3780° to 1740° R. 
However, if surface temperatures below 1000° R are required, radiation 
reduces the required convection heat transfer only slightly. 


CONCLUDING REMARKS 


Relations are presented for calciilating the forced- convection 
heat-transfer coefficient and local rate of heat transfer at both two- 
dimensional and axially symmetric stagnation points . These relations 
involve the local velocity gradient at the stagnation point and the 

dimensionless heat-transfer parameter Nu/V^^- 


An examination of various exact solutions shows that cooling a 
suirface at constant temperature without transpiration cooling can re- 
duce the parameter Nu//^Re^ by less than 10 percent from that of the 


insiilated surface. 


On the other hand, flow into the boundary layer through a porous 
surface tends to cause a large decrease in heat transferred to the 
svurface . This decrease is more pronounced as the siurface temperature 
is reduced. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 20, 1955 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 
A constant from Ug = Ax™ 

B constant from (t•^^ - tgp) = Bx^ 

C stagnation-point velocity gradient 

c specific heat at constant pressure 

ir 

f boundaiy-layer stream function 

(AH)^g enthalpy rise of transpiration coolant 


h 

k 

M 

m 


heat-transfer coefficient = 



thermal conductivity 
Mach number 

exponent from Ug = Axi™ 


Nu 

n 

Pr 

<1 

*3-ic 

Ir 

Re^ 

r 


Nusselt nixmber = — 

exponent from (t^ - t^p) = Bx^ 

Prandtl number 

local heat-transfer rate per unit area 
rate of heat transfer to an internal coolant 
radiation heat-transfer rate per iinlt area 

UgX 

Reynolds nimiber Re^ = 

'’w 

radius of sphere or cylinder 


I 
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S enthalpy function = r 1 

■'^0 

t temperature 

t^ vaporization temperatvire 

u velocity component in x-direction 

V velocity component in y-direction 

X coordinate from stagnation point along body surface 

y coordinate normal to body surface 

2m 

P press'ure- gradient parameter, p = 

X ratio of specific heats 

Tj boundary-layer similarity parameter 

p dynamic viscosity 

V kinematic viscosity 

p density 

Subscripts : 

aw adiabatic wall 

e local flow outside boundary layer (external) 

s stagnation value 

sp stagnation-point value 

w wall or surface value 

0 free -stream stagnation value 

1 conditions in front of body (but behind shock wave in super- 

sonic flow) 


r 
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APPENDIX B 


EXACT SOLUTIONS FOR STAGNATION-POINT FLOW WITH 
PRANDTL NUMBER OF 0.7 

In reference 9 Stewartson's transformation was applied to the 
compressible laminar -boundary- layer equations. The requirement of 
similarity was introduced and a viscosity law similar to that of 
Chapman and Rubesin (ref. 10 ) was ass\mied. The eqiiatlons that re- 
sulted are 


S" + Pr fS' = (1 - Pr) 


fi.. ^ pft. _ _ 1 _ s) 

(r - i)m2 


1 + 


r - 1 




(Bl) 


(f'f" + f"2) 


where f is the bo\mdary-layer stream function, S is an enthalpy 
function, and primes denote differentiation with respect to the 
boundary-layer similarity variable tj . For stagnation-point flow. 
Me = 0. With Pr = 0.7, the equations to be solved are 

f." ^ ^ p^j.,2 _ 1 _ s) 

S" + 0.7 fS' = 0 

with the boundary conditions 


(B2) 


f(0) =. 


^0 


without transpiration cooling 




w 


"•w 




vr for transpiration cooling 


f'(0) = 0 

S(0) = S^ = - 1 


> (B3) 


0 




For a two-dimensional stagnation point, 3=1; while for an axially 
symmetric stagnation point, 3 = l/2 (ref. 17 from Mangier 's trans- 
formation). Solutions to equations (B2) with boundary conditions 
(B3) were obtained on the IBM card- programmed calculator. 
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Values of the heat-transfer parameter Nu/^Re^ are obtained from 
the solutions as follows . For two-dimensional flow. 



so that, at a two-dimensional stagnation point (p = l): 

Nu ^ f 

For axially syrametric flow, the appropriate two-dimensional value 

of Nu/^Rey must be multiplied by 'y/s' according to Mangier ^s trans- 
formation, so that at an axially symmetric stagnation point (p = 1/2): 

The solutions of reference 7 are, of course, also adapted to axi- 
ally symmetric stagnation-point flow by obtaining the value of Ru//^Re^ 

for 3 = 1/2 (Euler number Eu ^ l/s) and multiplying by a/s". Since in 
reference 7 solutions are not presented for Eu = l/3, these must be 
obtained by interpolation. The normal velocity parameter defined in 

terms of the variables of reference 7 may be written 
^w/V^ = ■ (Eu + 1)/2'\/Eu, which for an axially symmetric stagna- 
tion point becomes v^//\/ = - 2W3. 
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(b) Axially symmetric. 


Figure 1. - Effect of Prandtl mamber and wall temperature on 
heat- transfer parameter. 


Heat- transfer parameter^ Nu/-^/Se^ 
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Normal velocity parameter, 


(a) Two-dimensional. 


(b) Axially symmetric. 


Figure 2. - Effect of transpiration cooling on heat-transfer parameter 
for various wall temperatures. Prandtl number, 0.7. 
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Convection or 'transpira'tion heat- Normal veloci ty p arameter^ 

transfer rate, q, Btu/(sec)(sq ft) v /VvC 

W' ^ w 

Figure 3. - Effects of internal cooling, transpiration cooling, and radiation on stagnation- 
point heat transfer and surface temperature of hemispherical-nosed body. 
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